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SUMMARY 

Yin  Chapter  I,  a  finite  element  methodology  is  described  for  analysis 
of  steady  nonlinear  periodic  oscillations  of  axisymmetric  shells.  The 
technique  is  a  combination  of  asymptotic  analysis  and  the  finite  element 
method  and  the  results  chat  it  yields  include  the  significant  effect  of 
nonlinear  interaction  between  the  waves  that  travel  in  opposite  directions 
along  the  circumference. 

In  Chapter  II,  an  extended  numerical  method  is  described  which  is 
applicable  to  shells  rotating  at  a  constant  speed  about  their  geometrical 
axes. 

In  the  final  Chapter  ILL*  a  general  solution  for  the  problem  of  non¬ 
linear  flutter  of  axisymmetric  shells  is  presented.  The  method  is  expected 
to  be  useful  for  developing  a  numerical  scheme  for  analyzing  this  class  of 
problems.  ,  /  ;  .  ■  -  rt  ''PAf.K  :fr)\  *  v l,  •  .<  >, 


CHAPTER  I 


FINITE  ELEMENT  ANALYSIS  OF  NONLINEAR  OSCILLATIONS 


OF  SHELLS  OF  REVOLUTION 


2 


INTRODUCTION 

In  this  chapter  a  computational  methodology  is  presented  for  analysis 
of  steady  nonlinear  harmonic  oscillations  of  stationary  shells  of  revolu¬ 
tion.  The  development  is  based  upon  the  results  which  were  derived  earlier  [1] 
in  the  context  of  the  more  general  problem  of  vibrations  of  slowly  rotating 
axisyrcnetric  shells.  In  brief,  the  results  in  [1]  are  in  the  form  of  a 
sequence  of  linear  equations  whose  solutions  can  be  used  for  the  calculation 
of  the  coefficients  in  the  asymptotically  exact  nonlinear  amplitude-frequency 
equations  for  the  shell.  In  what  follows  is  presented  a  finite  element  method 
for  solution  of  these  linear  problems;  also  described  herein  are  some  numerical 
results  which  have  been  obtained  from  an  implementation  of  the  proposed  scheme 
on  a  computer. 

In  the  class  of  problems  of  interest  in  this  chapter,  the  case  of  a 
circular  cylindrical  shell  has  been  most  extensively  studied.  A  historical 
account  of  various  investigations,  which  need  not  be  recited  here,  was  given 
by  Evensen  [2]  about  a  decade  ago.  Although  even  at  that  time  some  systematic 
perturbation  analyses  of  a  simply  supported  circular  cylindrical  shell  were 
available  [3,4],  further  studies  of  the  simpler  problem  of  a  circular 
ring  [1,5-7]  clarified  a  number  of  issues  arising  from  Evensen's  review.  It 
is  apparent  from  this  review,  and  from  the  results  in  [1,3,4]  that  the  most 
interesting  feature  of  nonlinear  ’scillations  of  shells  of  revolution  is  the 
interaction  between  the  circumferentially  traveling  waves  that  occurs  due  to 
the  rotational  symmetry  of  the  structure.  In  light  of  the  facts  that  some 
proposed  finite  element  based  schemes  for  analysis  of  nonlinear  oscillations 
of  elastic  structures  even  in  absence  of  mode-interaction  are  not  generally 
valid  [8]  and  that  attempts  to  use  the  finite  element  method  for  vibrations 
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of  axisymmetric  shells  (by  virtue  of  yielding  unacceptable  results)  have 
generated  inuch  controversy  [9-11],  the  methodology  presented  in  what  follows 
appears  to  be  of  some  technical  interest. 

Since  the  object  of  the  analysis  of  the  type  described  here  is  to  ob¬ 
tain  the  relationships  between  the  amplitude  and  frequency  of  oscillation 
of  the  shell,  the  same  calculations  could  be  done,  at  least  in  principle, 
by  using  a  commercially  available  finite-element  code.  Such  a  computation 
would  entail  a  series  of  time-integrations  with  different  set  of  initial 
conditions,  with  the  amplitude  and  frequency  obtained  from  the  time-historv 
of  notion.  This  methodology  is,  however,  expected  to  be  quite  uneconomical 
in  comparison  to  the  asymptotic  method  we  propose  here:  the  economy  evidently 
results  from  a  systematic  elimination  of  both  time  and  the  circumferential 
coordinate.  The  efficiency  of  the  scheme  becomes  more  evident  in  case  of 
the  forced  vibration  problems  for  which  the  response  of  the  shell  is  quali¬ 
tatively  very  rich  (cf.  Fig.  2  of  [3]).  Although  the  forced  response  curve 
on  the  basis  of  the  asymptotic  analysis  can  be  obtained  by  solving  a  small 
number  of  nonlinear  algebraic  equations,  it  is  unlikely  that  a  brute-force 
time  integration  scheme  can  predict  all  the  features  of  the  rather  compli¬ 
cated  behavior  of  the  shell  near  a  natural  frequency  as  in  [3]. 
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ASYMPTOTIC  RESULTS 

For  Che  sake  of  easy  reference,  in  this  section  we  present  the  results 
for  stationary  shells  obtained  from  [1],  For  our  purpose  these  results  have 
obviously  been  modified  to  exclude  the  terms  that  arise  due  to  the  rotational 
motion  of  the  shell. 

It  is  convenient  to  write  the  equations  of  motion  in  the  weak  form 

Mn(q,:q)  +  Pu(q,?q)  +  P21(q,6q)  +  P31(q,5q)  =  0  ,  (1) 

where  the  notation  due  to  Koiter  [12]  has  been  used.  Specifically,  the 
array  q  represents  the  three  displacement  components  of  the  shell.  The 
cuantitv  P  (q,5q)  in  (1)  is  the  first  variation  of  the  quadratic  term  in 
the  functional  for  strain  energy,  i.e. 

Pn(q,{q)  =  SP2(q)  (2) 

and,  similarly,  P^  and  P^  denote  quadratic  and  cubic  terms,  respectively, 
in  the  equilibrium  equations.  They  represent  the  first  variations  of  homo¬ 
geneous  functionals  of  degree  three  and  four,  respectively,  i.e., 

p2i(c*’5ci)  =  P31(q,5q)  =  <5PA(q)  •  (3) 

As  should  be  obvious,  the  first  term  in  (1)  represents  the  inertia  of  the 
shell  in  terms  of  the  first  variation  of  a  quadratic  functional  M^,  with 
M2(q)  being  the  kinetic  energy  of  the  shell. 

The  form  (1)  of  the  equilibrium  equation  represents  them  in  a  very  com¬ 
pact  fashion.  The  form  includes  both  the  differential  equations  and  the  natural 
boundary  conditions  (or  algebraic  equations  for  discrete  systems)  and  is 
valid  for  any  structure.  Moreover,  for  shells  of  revolution,  which  are  of 
interest  here,  each  of  the  functionals  in  (1)  is  of  the  type  exemplified  by 
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P,(q) 


2- 

f 

|  rP,(q)  d"  ds* 

r 

'  0 


(4) 


with  the  components  of  q  depending  upon  the  angular  coordinate  ”,  another 
snatial  coordinate  s,  and  time  t.  Evidently,  all  the  components  of  q  must 
be  periodic  functions  of  r ,  with  a  period  2". 

In  order  to  obtain  solutions  of  (1)  which  are  periodic  in  t:me  with  a 
frequency  ,  it  is  convenient  to  scale  rime  by  u>.  In  terms  of  the  scaled 
time  coordinate  the  equilibrium  equations  become 

c“Mn(q,<q)  +  Pn(q,5q)  +  P21(q,5q)  +  P^Cq.fq)  =  0,  (5) 

where,  now,  the  superposed  dots  denote  derivative  with  respect  to  the  non- 
dimensional  time.  To  (5),  one  must  append  the  appropriate  kinematic  con¬ 
ditions  on  the  two  ends  and  s2  of  the  shell,  together  with  the  periodicity 
conditions 


q  (s ,  9  ,  t )  =  q(s,6+2r,t)  =  q (s , 9 , t+2r ) .  (6) 

Equations  (5,  6)  define  a  set  of  nonlinear  equations  with  u  as  a 
parameter;  these  equations  have  a  trivial  solution  for  all  values  of  u>.  The 
condition  for  the  existence  of  bifurcation  points  on  this  trivial  solution 
branch  in  the  (u>,q)  space  leads  to  the  classical  linear  free  vibration 
problem 

-w02Mn(y.6q)  +  p11(y»<5c0  =  °*  (?) 


where  is  a  natural  frequency,  and  y  is  the  associated  natural  mode  of 
oscillation.  Due  to  rotational  symmetry,  there  are  two  linearly  dependent 
solutions  of  (7)  for  a  given  natural  frequency  u^,  i.e.,  y(s,e)  is  of  the 
form 


ine  .  -in© 

,  Y(s)e 


* 


Y  (s)e 


with  the  time  dependent  motion  being  a  linear  combination  of  the  four 


sc i ut ions 


.  .  tin-  tit 
\{s)e  e 


"his.  of  course,  is  the  result  based  upon  the  linear  theory.  The  asymptotic 


n-alvsie  in  [11  yields  that  near  a  natural  frequency  the  solution  of  (1) 


is  uiveo  by  (with  c.c.  denoting  the  complex  conjugate  of  the  term  preceding 


q  =  (jye’*'t  c.c.)  +  ( ■  ye”* 1  +  c.c.) 


+  oh„<:V“  +  cc.>+£-2  <0) 


c'))  o-it  _  rn) 

+  (urz^"  e“  +  c.c.)  +  (airZ-p  ’  +  c.c.)  +  h.o.t.. 


(S) 


where  h.o.t.  represents  terms  which  are  at  least  cubic  in  the  amplitudes  a  or 


(k)  . 


linear  in  U  -  ,g).  The  participating  modes  z„  in  (8)  are  obtained 


ij 


from  the  solution  of  the  linear  problems 


-M11(z11(2),5q)  +  P11(zu(2),fq)  +  P21(y,5q)  =  0 


11'  11 


(0) 


Pn  ( zi  i  >5ci)  +  pm(y,y,iq)  =  °» 


(9a) 

(9b) 


-Au02Mii(2i2(2)^q)  +  Pn(*12(2).*0  +  Pm(y.y^9)  =  0 


(0) 


pn<‘zi2  ,<5q)  +  piii^y»y«59)  =  o, 


(9c) 

(9d) 


with 


z  (2)  -  (2)  (0)  (0) 

Z22  _  Z11  *  22  11 


(10) 


(k) 


and  6-dependence  of  z. .  is  of  the  form 

ij 


(2) 


(2),  .  2in9  (0) 


'11  w,v/  “11  »  “n  '“’w/  “ii 


(».«)  -  i'*' z  „""<s,9>  -  Z„(0><.) 


(ID 


(2) 


(2) 


■n  (s,6)  -  Z  '‘'(s),  z.. 


(°)  _  7  (0),  V  2ine 

“i q  vS/C  • 
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Finally,  the  (complex)  amplitudes  a  and  £  are  related  to  the  frequency  u 
through  the  amplitude  frequency  equations 


2  2  2- 

-mO  -Wq  )i  +  a  +  y^ctrr-  =  0  , 

02a) 

2  2.  2  - 

-m(_  -*0  )  r  +  -,1::  r  +  =  0, 

(12b) 

where 


m  =  M11(y,y), 

(13a) 

'l  =  P111(Z11  ’  y,y)  +  P111('Z11  ),y,y)  +  P211(y,y’y)’ 

(13b) 

'2  =  ?lll(zll  ,y,y)  +  P111(Z12(  ),y,y) 

+?lll(z12(0),y,y)  +  2Poil(v’y>y)- 

(13c) 

Thus,  the  essential  effects  of  nonlinearities  on  vibration  of  an  axisymmetric 
shell  can  be  determined  by  calculating  the  coefficients  y  and  y^  in  (12) 
by  means  of  (13b, c).  As  a  result,  in  addition  to  the  solution  of  the  linear 
free  oscillations  problem,  it  is  sufficient  to  obtain  the  solution  of  the  four 
linear  equations  (9).  In  the  next  section  a  finite  element  method 
is  described  which  yields  discrete  approximations  to  equations  (7),  (9)  and 
(13),  and  this  approximation  allows  one  to  implement  the  scheme  just  presented 
on  a  digital  computer. 
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DISCRETIZATION 

The  development  in  this  section  is  based  upon  the  shell  theory 
derived  by  Sanders  [13]  under  the  classical  Kirchoff's  hypothesis  and 
the  small-strains,  moderate  rotations  approximation. 

We  first  list  the  explicit  form  of  various  functionals  in  (1).  If 
the  shell  is  assumed  to  be  made  of  an  orthotropic  material,  with  the 
principal  directions  of  orthotropv  being  along  the  circumferential  (•-) 

and  meridional  (s)  coordinates,  the  strain  energy  is  given  by  (see  Fig.  I-' 
'  s2  l2' 

P(q)  =  !  j  r[p2(q)  +  P3(q)  +  P4(q)]d0ds 
s  1  -  c 

-  P,(q)  +  P 3 (q )  +  P4(q)  ,  (14) 

where 


p2(q)  =  \  (Cues2s  +  C22£2e  +2C12€sseeS  +  C66EL> 


12  ss  96  66  str 


+  2  <Dn%s  +  b22<29  +  +  D66'se> 


(15a) 


■  Cll%s(Ss  +  ^  +  C22E99(Se  +  e2> 


+  C12|ESS<69  +  B2)+E69<6S2  +  B2)1 


+  C66Es9BsB6  ’ 


(15b) 


P4(q)  =  |  [Cn(sf  +  B2)2  +  C22(S2  +  82)2  +  2C12(bJ+62)(62+62) 


+  c66esee  3  • 


(15c) 


The  linear  strain  components  e  ,  eoq,  e  tc  ...etc.,  together  with  the 

SS  00  os  ss 

three  measures  of  rotational  displacements  8,8  and  6  are  given  in  terms  of 

s  6 
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the  three  translational  displacement  components  of  the  shell;  the  relevant 

equations  are  given  in  the  apnendix.  The  quantities  C. .  and  D. .  in  (15) 

ij 

are  evidently  dependent  upon  the  elastic  constants  of  the  shell  material 
and  the  shell  thickness. 

The  Linear  Free  Vibration  Problem 

The  discrete  approximation  of  the  linear  free  vibration  problem  (7) 
can  he  obtained  in  a  straightforward  manner.  The  displacements  associated 
with  the  natural  mode  of  oscillation  v  are  written  as 

>’  =  (U(n)(s),  -iV(n)(s)  ,  W(n)(s))ein',  (16) 

so  that  on  using  the  kinematic  relations  one  can  obtain  the  other  kinematic 
variables  associated  with  v  in  the  form 


(e  ,  e.  ,  e  „)  =  (E(n),  -i  E(^)elnS  , 

ss’  St’  St  SS  00  ’  s6 


(17a) 


(S  ,  3  ,  S)  =  (B^n)  ,  -i  B^n),-  i  B(n))einB  , 

S  S  u 


(17b) 


(<  <  <  )  =  (<(n)  ,c(n)  -  i  <(n)  )ein9 

Kss’  66’  <s9;  ^  ss  ’  *99  ’  1  s6  )e 


(17c) 


The  relationship  between  the  quantities  E^"^ ,  E^”\ . .  etc.,  and  displace¬ 
ment  variables  and  has  been  given  in  the  appendix. 

In  a  manner  similar  to  (16),  the  virtual  displacements  6q  in  (7) 
are  written  as 

CO 

(6u,  5v,  6w)  =  l  {(6U(n),  -i6V(n),  6W(n))ein6  +  c.c.}  (18) 

n*0 


with  analogous  expressions  for  the  associated  6E  ,  t$E._,  6E  _  ,  ...  etc. 

SS  00  so 

Substitution  of  these  expressions  into  (7)  with  (14,  15a)  yields  ^ 

■S2 

Pn(y,6q)  =  2tt  r[5£(n)T$  +  6<M\M]ds,  (19) 

_ S1 

^  The  symmetric  matrix  C  in  (19)  has  the  components  C  ,  i,j  *  2,  with  its 
third  column  being  (0,  0,  C66)T,  cf.  (15).  The  matrix  £  is  defined 
similarly. 


where  =  (E^“  ,  E"w  ,  )  ,  and  .similarly, <EV"' ,  *< uw  and 

are  arravs  with  real  elements.  One  can  now  choose  appropriate  one¬ 
dimensional  basis  Sanctions  for  the  displacements  =  (U^n\  V^n\  \’^n^ 

(cf.  (161),  i.e.  , 


.. (n )  /  \  _  pT  (n) 
JL  (s)  =  F  (s)  y 


(20) 


Cn ) 

where  f  represents  the  shape  functions  and  y  is  the  array  cf  generalized 


displacements.  rinallv,  with  the  relations  (A. 10,  12)  written  in  the  more 
com.nact  form 

E(n)  =  F(n)lJ(n),  <(n)  =  ?(n)U(n),  (21) 

and  with  (20),  we  obtain 

Pn(y,:<,)  :  :-,wV”>,w  , 


(22) 


where 


s(n)  -  j  [(|(nV)Tf/lfnV) 


(n)_TsT_ (n)_T.  . , 

+  <$  l  >  l  ^ds* 


(23) 


is  the  stiffness  matrix  associated  with  a  mode  of  deformation  with  the 
circumferential  wave  number  n.  The  corresponding  mass  matrix,  being  ob¬ 
tainable  from  the  functional 
So  r 2tt 


'^(y,  <5q)  = 


fs2  r 


ro(u6u  +  v£v  +  w£w)d9ds, 


(24) 


51  0 


can  be  calculated  in  a  similar  fashion,  and  is  given  by 


(n) 


rs  i 


rP^Tds, 


(25) 


where  p  is  the  density  of  the  material  of  the  shell.  Obviously,  both  of 
these  matrices  are  obtained  by  means  of  the  usual  finite  element  assembly 


process.  The  solution  of  the  generalized  eigenvalue  problem  associated 
with  these  matrices,  i.e.. 


<--'Vn>  +  K<n,)5-<n>  -  0.  ;!/">!«  4  0, 


(26) 


provides  an  approximation  to  the  solution  of  the  problem  (7),  the  linear 
free  vibration  problem. 


The-  Froblem.s  for  the  Par ticipating  Modes  and  the  Ampl itude-Frequencv  Equations » 
For  the  solution  of  the  nonlinear  problem,  it  now  remains  to  obtain 
discrete  analogs  of  (9a-d)  and  (13).  For  this  purpose  it  is  convenient 
to  first  consider  the  forcing  terms  in  (9)  (i.e.,  P  (v , 5q) »P-q^ (y ,y , - q) ,etc. ) 
which  require  the  calculation  of  the  first  and  second  variations  (Frechet 
derivatives)  of  the  functional  P^(q)  defined  by  (14)  and  (15b).  Although 
the  algebra  is  tedious,  the  calculations  are  conceptually  quite  straight¬ 
forward.  We  begin  by  noting,  with  equations  (3),  (14)  and  (15b),  that 


p21(q,<5q)  * 


•S2-2-1 

J  rp21(q,fq)dfds. 


and  that 


2„21(,.5q)  -  Jc„rcn(^  +  S2)  ♦  C12(B92  *  62)] 

*  <E69,C12(e82  *  e2)  +  C22<66  +  62)1 


+  5c  [C,,£  6  ] 
st>  66  s  6 


+  6£S[2C11CSS6S  +2ci2Eee6s  +  c66EseS9] 

+  f£9l2C12Css6e  +  2C22C99e9  +  C66ese6sJ 
+  66[2cii£sss  +  2C22EeeB  +  2c12(ess  +  cee)6l< 


(27) 


C28) 
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In  order  to  evaluate  P0^(v,6q)  for  (9a),  we  use  the  equations  (16,  17a, b, 


L 8 ,  27)  and  (28).  The  result  is 


P?1(y,-q)  =  2-  |  r(.'E(2n)Ta  +  5p(2n)TJ?)ds, 


where 


31  =  [Cll(Bin)  “  B(n)  ]  _C12(B-n)  +  B(n)  )]/2’ 


(30a) 


a^  =  [C  (B(n)  -  B (n)  )  -  C  (B,(n)  +  B(n)  )]/2, 

2  s  22  5 


(30b) 


*3  '  C66 


(30c) 


b  =  C  E(n)B(n)  +  C  E^Vn)  -  C,,E(n)B!:n)/2, 
1  11  SS  S  1  "  “fie  ct  S 


b2  -  C12tss)Be">  +  +  C^)B<")/2, 


22  55  9 


66  s5  s 


(30d) 


(30e) 


b3  ’  [CllEs^  +  C22E'”)  +  C12(Eis>  +  E‘")”B<“)  ' 


( 30f  ) 


Ve  recall  here  that  the  quantities  etc.  in  (30)  are  related  to 

the  displacements  associated  withy  via  equation  (A. 10,  11). 

The  discrete  analog  of  (29)  is  obtained  by  using  (18,20,21).  Thus, 
first  the  quantities  a^  and  b^  in  (30)  are  calculated  from  (A. 10,  11)  and  (18). 
For  example , 

a3  =  CW^,V  0:  ~ H  [0;  sini>/r:  n/r]£T(s)y )  . 


From  the  solution  of  the  homogeneous  problem,  the  quantity  a^  and, 

similarly,  all  components  of  ^  and  £  can  be  calculated  for  any  value  of  s. 

Using  these  functions  in  (29),  one  obtains 

T 

P21(y,6q)  2:  2*  5q(2n)  £(2)  ,  (31) 


where 
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,'s2 

,(2)  i  ,  ‘(2n)_T,  ,,T  ,  ,  .  r_(2n)_T,  ,.T. 
f  =  r- IE  F  (s' ]  $(s)  +  [g  F  (s)J  Jt>(s)}  ds 

I  ^  c 

'S1 


(32) 


In  ( 3 2 >  we  have  used  the  definition  (21)^  and  an  analogous  definition 


B(l)  =  B(l\  B(1))T  =  B(l)U(l), 

x  s  -  x  x  ’ 


(33) 


2(1) 


with  B  obtained  from  the  equations  (A. 11). 


The  remaining  forcing  term  in  (9)  is  P^^(v,y,;q)  for  which  the  second 
variation  of  needs  to  be  calculated.  Since  the  details  are  similar  to 
the  ones  described  for  computing  we  list  only  the  final  results. 

Thus,  one  obtains 

rs2 

.T  ,  J 


Pin(y,y,iq)  =  jg  r(:£(o)  c  +  5B(o)  £)ds. 


where 


C1  =  C11  <B^  +  B(n)  }  +  C12(B0n)  +  s(n) 


c  =  C  (B(n)  +  B(n)  }  + 

l  xz  s 


c22(Sgn)  +  B(n)  ), 


(34) 

(35a) 

(35b) 


c3  =  0, 

d  =  2  C  E(n)B(n)  +  2C.,E^)B(n)  +  C  E("Vn)  , 
1  11  ss  s  12  60  s  66  s6  0 


(35c) 

(35d) 


d2  =  d3  =  0. 


(35e) 


Finally,  the  discrete  analog  of  (34)  turns  out  to  be 
"ill1 


pm(y.y»*q)  1  2tt 


with 


l 


(o) 


r {  [£(CV(s)]T  £  +  [|(0)FT(s)]^}ds. 


(36) 


(37) 


Our  description  of  the  procedure  for  discretization  of  (9)  is  now 
essentially  complete.  With  the  vectors  and  defined  by  (32)  and 


(37),  the  approximations  to  the  participating  modes  z„  in  (9)  can  be  ob¬ 
tained  by  solving  the  linear  systems ^ 


<-s¥2n)  +  ¥(2n))«g)  ♦  /(2) 


-  0, 


(o)„  (o)  (o) 

K  *11  +  1  “  °’ 


(38a) 

(38b) 


/  ,  2m(0)  ,  ..<0)._(2)  ^  f(0) 

(--.  o  M  +  4,  )Z12  +  f  =  0  , 


(38c) 


(2nW0)  ,  or(2) 


r~"'i HY  +  2V  =  °- 


(38d) 


In  (38)  M^°\  anc*  are  defined  by  (25,23)  with  n  replaced 


(k)  (k) 

by  zero  and  2n,  respectively,  and  Z_  are  discrete  analogs  of  Z  in  (9). 


Once  the  participating  modes  have  been  calculated  from  (38),  one  only 
needs  to  use  (13)  to  calculate  the  coefficients  in  the  amplitude-frequency 
equations  (12).  All  the  terms  in  (13)  except  P2^^(y»y>y)  can  be  obtained 
by  using  equations  (31)  and  (37).  Consequently,  we  proceed  to  calculate 
the  last  term  in  (13b).  For  this  purpose  we  write  (15c)  in  the  form 


p4(<i)  -  llCnOjV  +  c22qTy2  +  2c12(,Ta)<rt> 


+  c66eJeel- 


(39) 


where 


r,.(6s:8)T,  £  -  (B,:  S)T. 


(40) 


The  second  variation  of  P^(q)  is  then  found  to  be 


I 


We  note  that  by  definition. 


pni(y.y.<5q) 


2P2i(y,5q)  so  that  the  forcing 


term  in  (9d)  is  twice  that  in  (9a). 
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r s-,  , 2- 

!  “  I 

P^jCq,  q^'q)  =  j  !  rp-.-^q,  q^qldsde 
-sr°  - 


where 


2p211fq,  qlt:q>  =  ^  x  [  2  ( 'T' .  )  (:7  ?;.)  +  (nT:.)  (:  ^p)  ] 

+  c22i2((Ti1)aTH> +  (fVujipi 

+  Cpi’fJ^r^n)  +  (;TU(-T- ,) 

+2(n  1r)(^T?)  +  (rTn)({CTC  )] 

->  Av  %  AV  '"'w  'V  '"w  '  v 

+  4  C*J<=  +  5s-£  ,)i  S- 

2  66  -Is  -  si  s  ? 


(41) 


(••■..-s  +  ':?.-fs)(V61  +  fs1  =  0)1* 


(42) 


and  the  kinematic  quantities  with  subscript  1  are  associated  with  displace¬ 
ments  corresponding  to  q^.  We  now  use  (17b)  and  (40-42)  to  obtain 


P211(y,y,y)  =  2 ^ 


p  r  ds  =  2~y 


v’here 


p  =  3tC11(2g12  +  &22)  +  C22(2h12  +  h22)] 

+  C12(2glhl  +  82h2)  +2  C66^B0  Bs  ^ 


and 


,  -  Bs^>2  +  .<<  g2  -  Bs^>2  -  B(n>2 


(43) 


(44) 


(44a) 


hl  -  (B0(n)  +B(n)2),  h2  =  -  hj.  (44b) 

Thus  once  the  free  vibration  problem  has  been  solved,  the  left  hand  side  of  (43) 
can  be  evaluated  by  using  the  equations  (20),  (A. 11)  and  (43).  Finally,  from 
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In  summary,  the  computations  proceed  in  the  following  manner.  First, 
for  a  given  circumferential  wave  number  n,  the  stiffness  and  mass  matrices  are 
obtained  by  using  (23,25).  The  solution  of  the  generalized  eigenvalue 
oroble..<  (26)  yields  a  natural  mode  of  oscillation  and  the  associated 
frequency.  The  natural  mode  is  used  for  calculation  of  the  forcing  terms 
f^  and  via  (32,37).  Then  the  mass  and  stiffness  matrices  corres¬ 

ponding  to  circumferential  wave  numbers  2n  and  zero  are  assembled  so  that 
the  participating  modes  can  be  calculated  from  (38).  Finally  with  the  scalar 
y  calculated  by  means  of  (43,44),  equations  (45a-c)  are  used  to  compute 
the  desired  coefficients  the  amplitude-frequency  equations  (12).  The 
results  from  a  computer  implementation  of  the  procedure  just  described  are 
presented  in  the  next  section. 


NUMERICAL  RESULTS 

In  the  program  that  has  been  developed  on  the  IBM  4341  at  the  Yale 
Computer  Center,  one  dimensional  Hermite  cubics  were  used  as  the  shape  functions 
for  all  the  three  displacement  components,  together  with  four  point 
Caussian  quadrature  rule  for  numerical  integration.  All  the  computations 
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were  done  in  doable  precision.  For  the  purpose  of  comparing  the  results 
with  those  obtained  previously,  the  scalars 

?1  =  /(8mw  2),  (46a) 

P0  =  (,2  -  Yj)  '(16hl.-  2 )  ,  (46b) 

were  computed.  The  scalars  occur  in  a  set  of  amplitude-frequency  equations 


equivalent  to  (12),  to  wi t , 

-(22  -  l)a  +  (P  +  P2)a3  +  (P  +  P2cos2-Dab2  =  0,  (47a) 

2 

ab  sin2v  P2  =  0,  (47b) 

-(f.2  -  l)b  +  (?1  +  P2)b3  +  (P1  +  P2cos2Dba2  -  0,  (48a) 

-ba2sin2(  P_  =  0  ,  8  =  w/w  .  (48b) 

l  o 

These  equations  are  obtained  by  setting 

a  =(a  +  ibe^v)/4,  (49a) 

£  -  (a  -  ibeili')/4,  (49b) 


in  (12),  adding  (12a)  to  (12b)  and  subtracting  (12b)  from  (12a),  using  the 
definitions  (46),  and  satisfying  both  the  real  and  the  imaginary  parts  of 
the  resulting  equations.  Evidently,  equations  (48)  are  equivalent  to  those 
derived  by  Ginsberg  [3]  for  the  special  case  of  a  circular  cylindrical  shell. 

In  Table  1  are  presented  the  results  of  a  convergence  study  for  a  simply 
supported  circular  cylindrical  shell  analyzed  in  [3].  Only  half  of  the 
shell  was  analyzed,  with  symmetry  conditions  imposed  at  s  *  0,  (w1  ■  v'  *  u  *  0) 
and  the  simple  support  conditions  (v  •  w  ■  0)  at  s/R  *  L/2R,  with  L  being  the 
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length  of  the  shell  and  R  its  radius.  It  is  evident  from  Table  1  that 

the  convergence  of  the  various  coefficients  with  increase  in  the  number 
of  elements  is  rapid  and  monotonic. 

In  Table  2,  the  computed  results  have  been  compared  with  those 
in  [3].  Although  one  of  the  reasons  for  the  differences  could  be  the 
different  shell  theorv  used  in  [3],  the  rather  substantial  discrepance 
cannot  be  explained  solely  on  this  basis.  (We  note  here  that  ultimately 
Ginsberg's  analysis  is  also  a  numerical  scheme  involving  summation  of  a 
convergent  series).  It  may  also  be  noted  here  that  a  direct  numerical 
counterpart  of  Ginsberg's  analysis  would  entail  the  rather  uneconomical 
computation  of  a  large  number  (theoretically  all)  of  the  natural  modes  of 
oscillation  and  the  associated  frequencies. 

In  Table  3  we  have  presented  some  results  obtained  for  a  simply 
supported  frustum  (conical  shell).  Specifically,  only  the  circumferential 
and  transverse  displacements  of  the  shell  were  constrained  to  vanish  at  its 
two  ends.  For  the  semi-vertex  angle  a  =  tt/2,  the  frustum  degenerates  to 
an  annular  plate,  and  our  results  indicate  that  for  this  case  the  non¬ 
linearity  is  of  the  hardening  type  for  both  the  circumferentially  traveling 
waves  and  standing  waves. ^ 

^  it  may  be  noted  here  that  from  (47-48)  the  following  amplitude  frequency 
equations  can  be  derived: 

(a)  b  =  0,  (Q2  -  1)  =  (P1  +  P2)a2;  Standing  Waves, 

(b)  a  =  b,  \p  «  0,  (ft2  -  1)  =  ZP^a2;  Traveling  Waves. 

The  character  of  these  solutions  can  be  readily  verified  by  using  (49)  in 
the  leading  part  of  (8). 
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CONCLUDING  REMARKS 

A  numerical  method  based  on  asymptotic  results  has  been  developed 
for  analysis  of  steady  nonlinear  harmonic  oscillations  of  shells  of  re¬ 
volution.  Some  results  have  been  described  for  circular  cylindrical 
and  conical  shells.  As  discussed  in  [1],  once  the  coefficients  in  the 
amplitude  frequency  equations  have  been  obtained  for  the  free  vibration 
problem,  the  problem  of  forced  vibration  can  be  reduced  to  solution  of 
two  complex  nonlinear  equations.  Thus  our  results  can  be  directly  used  for 
the  calculation  of  the  response  of  harmonically  excited  shells  of  revolution. 
Finally,  further  extension  of  the  numerical  method  to  include  the  effect  of 
steady  rotation  of  the  shell  is  the  subject  of  the  next  chapter. 
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APPENDIX 

A.  Strain-Displacement  and  Rotation-Translational  Displacement  Relations 
(See  Fig.  1  for  definitions  of  r,  r.  and  f). 

Membrane  Strains 


u  +  w  /  r 

(A.  1) 

>s  ; 

(v  /r  +  u  cos : /r  +  w  siny/r) 

>  - 

(A.  2) 

(v  +  u  a/r  -  v  cos^/r)/2 

5  S  j 

(A.  3) 

Rotations 


S  =  -  w  +  u/r 
S  ,s  $ 


=  (— w  .  +  v  sini)/r 


S  =  (v  +  v  cosi/r  -  u  /r)/2 

,  S  ,  D 


Bending  Strains 


<  =  -  w  +  u  /r.  -  u  r,  /r. 

SS  ,SS  ,S  <t>  (p,s  <1 


K09  =  (6e,e  +  6S  cos<p)/r 


:s9  =  lB9,s  +  Ss,9/r  '  E9  cos*/r  +  (sin*/r  -  i;VB]/2 


(A.  A) 

(A.  5) 
(A.  6) 

(A.  7) 
(A.  8) 
(A.  9) 


B.  Kinematic  Relations  for  the  Quantities  in  (17)  for  Deformation  Associated 


with  Circumferential  Wave  Number  n 


E(n>  ^ 
ss 

l  d/ds 

0 

> 

1/r 

♦ 

fD(n> ^ 

F(n) 

E00 

=* 

cos4>/r 

n/r 

siniji/r 

VM 

\  Ese  , 

%  * 

^-n/2r 

-|(d/ds-cos$>/r) 

0  1 

*"Vn) 


(A. 10) 


*-  V 
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Table  1 


Convergence  Studv  for  a  Circular  Cylindrical  Shell 
of  Isotropic  Material,  Poisson's  Ratio  v  =.3. 


No .  of 

I./R(1) 

n(2)  | 

1 

Elements 

J3) 

"o 

w(4) 

max 

p(5) 

V2  1 

5 

.6531 

.  9576 

-5. 361 

-74.17  j 

1 

2 

!0 

.  6531 

.9570 

-5.406 

-74.21 

1 

20 

.  6531 

.  9570 

-5.410 

-74.22 

5 

.2436 

.9914 

-565. 6 

-1271.3 

8 

10 

.  2436 

.9914 

-566.  7 

-1271.4 

1 

-  - 

20 

_ 

.2436 

.9914 

-566.7 

-1271.4 

5 

.1130 

.8836 

-1.489 

-61.70 

2 

.1130 

.8836 

-1.490 

-61.70 

4 

20 

.1130 

‘  ' 

.8836 

-1.490 

-61.70 

5 

.1824 

.9921 

-6.805 

-822.4 

8 

10 

.1824 

.9921 

-6.809 

-822.5 

20 

.1824 

.9921 

-6.810 

-822.5 

1.  L/R  =  Length  to  Radius  Ratio. 

2.  n:  Circumferential  Wave  Number. 

2  2  1/2- 

3.  a  *  [ (1-v  )pR  / E]  m  where  m  is  the  lowest  natural  frequency  associated 

o  o  o 

with  n  and  E  is  the  Young's  Modulus. 

=  w  /R,  where  w  is  the  maximum  normal  displacement  associated 
max  max  r 

2  2  2  ^ 

with  the  normal  mode,  normalized  such  that  (w  +v  +  u  )/R  ■  1, 

max  max  max 

as  in  [3]. 


5.  Coefficients  in  amplitude  frequency  equations  (47,  48). 


Table  2 

A  Comparison  of  Numerical  Results  with  those  in  [3]. 
(Parameters  and  Notations  are  same  as  in  Table  1.) 


New 

Results 


Results 
from  [3] 


New 

Results 


Results 
from  [3] 


New 

Results 


Results 
from  [3] 


New 

Results 


Results 
from  [3] 


"o 

w 

max 

■■ 

.  6531 

.  9576 

.6531 

.9576 

.2436 

.9914 

.2437 

.9914 

.1130 

.8836 

.1130 

.8836 

.1824 

.9921 

.1824 

.9922 

•74.22 


-96.00 


1271.4 


-1771 


-61.70 


-68.15 


.9921  -6.810  -822.5 


-873.1 


Table  3 


Coefficients  in  the  Amplitude  Frequency  Equation  (47) 
for  a  Simply  Supported  Frustum. 


- r 

l 

(Degrees) 

P1 

P2 

15  | 

i 

0.5991 

-6.49 

- - - - 

-57.71 

30 

0.  4983 

+1.  70 

-37.24  : 

45 

0. 3753 

2.85 

133.85 

60 

0.2511 

10.43 

-97.08 

75 

0.1303 

32.56 

-311.16 

90 

0.0340 

469.14 

233.60 

(Annular  Plate) 

[v  *  0.3,  R/h  =  100,  L/R  =  1,  (see  Fig.  2),  n  =  2.  Notation  is  the  same  as  in 
Table  1. ] 
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INTRODUCTION 

In  this  chapter  we  describe  a  finite-element  methodology  for  analysis 
of  steady  nonlinear  harmonic  oscillations  of  rotating  axisymmetric  shells.  As  in 
Cyan ter  1,  the  methodology  is  based  on  the  asvmptotic  results  derived 
earlier  [1],  and  it  yields  an  asymptotically  exact  nonlinear  relationship 
between  the  amplitude  of  oscillation  of  the  shell  and  its  frequency. 

The  main  influence  of  rotation  on  axisymmetric  shells  is  that  it 
destroys  the  reflection  svmnetrv  with  resoect  to  the  circumferential  direction. 

As  a  result,  the  two  waves  that  travel  in  opposite  directions  along  the 
circumference  are  expected  to  exhibit  different  behavior,  both  in  the  linear 
and  in  the  nonlinear  range.  It  is  for  the  assessment  of  such  differences 
between  the  two  circumferentially  traveling  waves  that  the  numerical  method 
developed  in  this  chapter  can  be  utilized. 

Although  in  mechanics  literature  there  is  a  large  number  of  accounts 
of  analytical,  numerical  and  experimental  studies  on  nonlinear  oscillations 
of  stationary  shells,  the  author's  analysis  [1]  appears  to  be  the  only  one 
in  which  explicit  results  have  been  obtained  on  the  oscillations  of  a  rotating 
shell  type  of  structure.  Due  to  the  analytical  intractability  of  the  general 
problem,  the  study  in  [1]  was  restricted  to  the  simple  problem  of  a  ring: 
evidently  the  finite-element  method  described  here  can  be  utilized  for  more 


general  problems. 
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r  A  RELATION 

As  in  Chapter  I,  we  use  the  Sanders'  strain-displacement  relations 

under  the  snail-strains,  moderate  rotations  approximation.  The  strain  energy 

(T ) 

of  the  shell  is  then  given  bv  the  functional  p(q)  defined  by  (1.14,  15)  *  . 

In  terns  of  this  functional,  the  equilibrium  equation  of  the  shell  in  a  co¬ 
ordinate  f rane  rotating  with  the  shell  can  be  written  as 


i?(q)  + 


c  (a  iu 
s 


+  a:v  +  a  iw)rdrds. 
r 


(1) 


In  this  equation,  n  is  the  mass  per  unit  area  of  the  shell  and  the  quantities 

...  2  2 
a  =  u  -  2.'v  cos:  -  T:  cosi(u  cos;  +  w  sini)  -Or  cosi  ,  (2a) 

s 

a.  =  v  +  2 Du  cos;  +  2Dw  sinj  -  fl^v,  (2b) 

••  .  2  2 
af  =  w  -  2?v  sin;  -  D  sini(u  cos<J>  +  w  sin$)  -Hr  sirni  (2c) 

are  the  acceleration  components  of  the  shell  in  the  rotating  frame,  see 
Figure  1.  Evidentlv,  the  quantity  fi  in  (2)  is  the  speed  of  steady  rotation 
of  the  shell  about  its  geometrical  axis. 

The  last  terms  of  the  equations  (2a, c)  lead  to  an  effective  centrifugal 
force  on  the  shell,  so  that  prior  to  analysis  of  oscillations  one  has  to 
calculate  the  steady  prestressed  state  about  which  the  vibrations  occur. 

This  prestressing  due  to  steady  rotation,  denoted  by  q°,  is  governed  by  the 
nonlinear  equations  obtained  by  dropping  from  (1,2)  the  terms  that  contain 
time  derivatives  of  displacements;  to  wit: 


(I) 


In  this  chapter,  equation  numbers  beginning  with  Roman  I  refer  to  the 
corresponding  equations  in  Chapter  I. 
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nn(nc,;c)  +  ?}  F11(f,fp1)  =  o. 


where 


2  f- 


Rn(q,  q)  = 


cos;(u  cos;  +  w  sin;):u 
+  v:v  +  sin*(u  cos;  +  w  sin:)‘w 


rdSds,  (Aa) 


r  ( f  r  )  =  _ 

n  ’  ■ 


(r  cos;  5u  +  r  sin;  fw)rd"ds. 


The  forcing  term  (Ah')  in  the  equation  for  prestress  due  to  rotation  is  ob¬ 
viously  ■— independent :  the  prestress  state  is,  therefore,  axisynmetric ,  as 
exnected. 

ror  a  fixed  value  of  T,  the  solution  of  discrete  analog  of  (3)  can  be 
obtained  bv  the  usual  Newton-Raphson  method  for  solution  of  nonlinear 
algebraic  equations.  We  write  the  solution  of  (1)  in  the  form 

q  =  q°  +  q  (5) 

% 

with  q  representing  the  oscillatory  motion  of  the  shell  superimposed  on  the 
steady  deformation.  On  substitution  of  (5)  in  (1,2),  together  with  the  use 
of  the  fact  that  q°  satisfies  (3),  we  obtain  the  homogeneous  equation  for  the 
free  nonlinear  oscillations  of  the  shell.  This  equation  can  be  written  as 


M^q.Sq)  +  Gn(q,6q)  +  Vn(q,5q)  +  V21(q,6q)  +  V^Cq.Sq) 
where 


2  r  2  it 


M_  ,  (q,6q)  =  I  (u5u  +  v5v  +  w5w)rd6ds, 
11  JSlJo 


f2lT  f  + 

G  (q,5q)  -  2fl  )  l 

si  o  L  +^v "  ^ 


6v)cos$ 


rd6ds, 


v5w)sin<f 


» 
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v11(q,tq>  =  T>11(q",rq't  -*■  1’111(q°,q.^)  +  P011  (q°,q  ,  ;a)  +  R^Cq’.iq),  (7c) 


V,1(q,-q)  =  p  2  ( q  , -q)  +  11  o  L  ]  f  »  q°.5a): 


V31(q,:q)  =  P^  -  >  “q > •  ('7e') 

Here,  as  usual,  F^(q,.q)  denotes  the  first  variation  of  the  quadratic  part 
of  the  strain  energv  functional;  P?^  and  p  are  defined  similarly.  We  note 
here  that  the  bilinear  functional  in  (7b)  is  antisymmetric  with  respect 

to  the  argument,  i.e.,  G^(q^,  92^  =  ~  S.1^2’  ^1^’  3S  C3n  reac*a^  verified. 

For  subsequent  discussion  it  is  convenient  to  dron  the  over-tilda  from  q 
and  to  scale  time  by  using  the  (yet  unknown)  frequencv  w  of  periodic  vibration 
of  the  shell.  The  equation  (6)  is  then  modified  to 


*  Mn(q,5q)  +  w  C^Cq.dq)  +  Vn(q,5q)  +  V21(q,iq)  +  V31(q,Sq)  =  0,  (S) 

where,  now,  dot  represents  derivative  with  respect  to  the  scaled  time  coordinate. 
With  this  scaling,  the  object  of  our  analysis  reduces  to  obtaining  solution  of 
(8)  whose  time  dependence  is  periodic  with  a  period  2 tt . 

An  asymptotic  solution  of  the  problem  has  been  presented  in  [l  ].  Here 
we  summarize  the  main  results. 

First  one  calculates  a  natural  frequency  u>q  and  the  corresponding  natural 
mode  of  oscillation,  y,  from  the  quadratic  eigenvalue  problem 

~<*>o^ii  (y » ^q)  +  ia)0  Gn(y.5q>  +  V11(y,5q)  =  0  ,  (9) 

to  which  must  be  appended  an  appropriate  normalization  condition  for  y.  If 
for  a  given  natural  frequency  the  associated  mode  of  oscillation  is  unique 
up  to  a  scalar  multiplier,  the  results  in  [1]  indicate  that  (8)  admits  solutions 
of  the  form 


( 
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q  =  (ave^  +  c.c.)  +  (  u“z  ^e^1"  +  c.c.)  +  aa  +  9(32,ol(u—w  ))  ,  (10) 

o 

where  c.c.  represents  the  complex  conjugate  of  the  term  preceding  it. 

(2)  (0) 

In  (10),  the  "participating  modes"  z  and  z  are  obtained  from  the 
linear  time-independent  problems 

-S2Mn(z(2),:q)  +  2i .  q  G11(z(:),fq)  +  Vn(z(2),:q)  +  V^v.f.q)  =  0,  (11) 


So) 


vn(zw,fq)  +  Vnl(v  ,y,;q)  =  0 


The  comnlex  amplitude  a  and  the  frequency  y  are  related  by  the  nonlinear 
amplitude  frequency  equation 

-(•>■"  -  -o“)ma  +  -  w Q ) q a  +  ^  a" a  =  0,  (13) 

where  the  nonlinearity  coefficient  y^ ,  and  the  other  scalars  are  obtained  from 


m  =  Mn(y*  y)>  g  =  i  y)> 


>1  =  2P?1(y,z(2))  +  P111<y1y»z(°))  +  p211(y»^*y)' 


The  main  object  of  our  calculations,  of  course,  is  to  obtain  the 
quantities  m,  g  and  y^  in  (13)  so  that  at  least  an  asymptotic  relationship 
between  the  amplitude  and  frequency  is  explicitly  available.  For  this  purpose 
as  we  have  discussed  in  this  section,  one  needs  the  computation  of  (a)  the 
steady  prestress  due  to  rotation,  (b)  a  natural  frequency  and  the  associated 
mode  of  oscillation  and  (c)  the  participating  modes.  Following  section  is 


devoted  to  some  aspects  of  these  computations. 
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THE  DISCRETE  ANALOGS 

Once  the  variations  of  the  various  functionals  needed  in  (9,11,12,14,15) 
have  been  calculated,  standard  finite  element  techniques  can  be  used  to  obtain 
the  discrete  approximations  of  these  equations.  Since  the  details  are  similar 
to  those  given  in  Chapter  I,  here  we  only  outline  some  of  the  main  results. 

On  discretizations  the  equation  (3)  for  prestressing  due  to  steady 
rotation  becomes 


,  /.  (o) . 
n(q  ) 


(o)  (o)  ,  1  ,  (o)  (o)  1  (o).  (o) 

h  q  +  X  Ni  (q  )q  +  T  NT2(q  )q 

\  L  *  r  r,  j  U  n, 


+  :.y°Vo)  +  .;2f  =  o 


where  K'  is  the  constant  stiffness  matrix,  with  N1  and  "2  being  symmetric 
matrices  whose  elements  are  linearly  and  quadratically  dependent,  respectivelv , 
on  tne  argument  the  vector  of  generalized  nodal  displacements.  Equation 

v 

(16)  can  be  solved  by  the  Newton-Raphson  scheme  (i  =  1,2,...£): 


with  £  determined  by  means  of  a  suitable  convergence  test. 

After  the  prestressed  state  due  to  rotation  has  been  obtained,  one  can 
solve  the  discrete  analog  of  the  linear  free  vibration  problem  (9)  which  turns 
out  to  be  of  the  form 


» 


1 


36 


Vn)  +  v  G(n)  +  K(n)]y(n)  =  0, 

An.  rv  "•  4  * 


(18a) 


(n)T  (n) 

\  X,  = 


(18b) 


where  n  is  the  circumferential  wave  number,  and  the  symmetric  matrices 


M^n\  and  K^n'  represent  the  mass,  the  Coriolis  acceleration ^ ^  and 


the  stiffness,  respectively. 

Although  some  attention  has  been  given  in  the  literature  to  the  development 
of  special  numerical  methods  for  solving  quadratic  eigenvalue  problem  of 
the  type  (18),  we  use  a  simple  Newton-Raphson  scheme  for  (18)  as  well;  to 


i  =  1,2,...,  till  convergence, 


i,2  (n)  i_(n)  „(n) 

-(w  )  M  +  w  G  +  K 

o  ^  o  v  v 


(~2w  iM(n)  +  G(n))y(n)i 

o  m  m 


,(15a) 


A,'")  +  „  G<”>  +  K<"))y<n)i 

o  'V  \  % 


v(n)i  (n)i  _  1 


(19b) 


We  note  that  if  we  write  the  displacement  components  in  the  form  (1.16) 
and  use  the  definition  (7b),  the  functional  iG^  in  (9)  becomes  a  real 
symmetric,  bilinear  functional. 
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=  (2n)if(2n)  +  «Of(o)  + 


(21c) 


where  the  last  terra  is  calculated  as  in  (1.43,  44). 


y.'^-ICAL  RESULTS 

The  algorithm  described  in  the  previous  sections  has  been  implemented 
on  the  iRv  414]  at  the  vale  Computer  Center.  The  program  employs  one- 
dimensional  Herraite  cubics  for  shape  functions  for  all  the  three  displace¬ 
ment  components  together  with  four-point  Gaussian  quadrature  scheme  for 
integration. 

In  Table  1  are  presented  the  results  for  a  simply  supported  circular 
cvlindrical  shell.  All  the  commutations  were  done  by  using  ten  elements 
for  half-lenath  of  the  shell  with  symmetry  conditions  imposed  at  the  mid¬ 
length.  It  is  evident  that  due  to  rotation  there  is  further  "softening"  of 
the  nonlinear  response.  However,  it  appears  from  the  results  that  in  the 
range  of  rotational  speeds  considered,  the  effect  of  rotation  on  the  non- 
linearitv  coefficient  y^  is  similar  for  both  the  forward  traveling  wave  and 
the  backward  traveling  wave.  Similar  trends  can  be  observed  in  Table  2,  in 
which  the  results  of  a  similar  calculation  for  a  circular  shell  clamped  at  both 
ends  have  been  presented. 
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Table  1 


Coefficients  in  the  .Amplitude  Frequency  Equation  (13) 
for  a  Simply  Supported  Circular  Cylindrical  Shell 
(L/R  =  1,  v  =  0.3,  R/h  =  100,  n  =  8) 


1 

i 

r  ^ 

i 

I 

r"  .  ~ 

“ 

g 

” 

! 

1 

I  o 

1 

F,B 

0.  243 

0.250 

0 

-  67.1 

1 

0.05 

F 

0.245 

0.250 

-  0.76  x  10~3 

-  67.9 

0.05 

B 

0.248 

0.250 

+  0.76  x  10‘3 

-  67.9 

0.10 

F 

0.251 

0.248 

-  0.15  x  10~2 

-  71.3 

0.10 

B 

0.257 

0.248 

+  0.15  x  10-2 

-  71.4 

0.15 

F 

0.262 

0.247 

-  0.23  x  10-2 

-  78.4 

0.15 

B 

0.272 

0.247 

+  0.23  x  10-2 

-  78.6 

0.20 

F 

0.279 

0.244 

-  0.30  x  10*  2 

-  90.2 

0.20 

B 

0.290 

0.244 

+  0.30  x  10'2 

-  90.4 

^  All  the  quantities  are  nondimens ional  with  the  scaling  parameters  being 
the  same  as  in  Table  1  of  Chapter  1. 


The  quantity  u)g  is  the  dimensional  lowest  natural  frequency  for  circum¬ 
ferential  wave  number  n  ■  8, 

F  denotes  forward  traveling  wave  (in  the  same  direction  as  the  sense  of 
rotation  of  the  shell)  and  B  denotes  backward  traveling  wave. 
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CHAPTER  III 

NONLINEAR  FLUTTER  OF  SHELLS  OF  REVOLUTION 


( 
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INTRODUCTION 

In  this  chapter  we  present  an  analysis  of  nonlinear  flutter  of  shells 
of  revolution.  A  comprehensive  summary  of  the  related  literature  can  be 
found  in  the  excellent  reviews  by  Dowell  [1,2]  and  Fung  [3]. 

A  number  of  experimental  studies  of  flutter  of  cylindrical  shells 
were  conducted  by  Olson  and  Fung  [4,5],  Stearman  [6]  and  Stearman  et.  al  [7]. 
Although  there  is  not  much  disagreement  on  the  qualitative  aspects  of  flutter 
between  the  observations  and  analytical  results,  the  agreement  on  the  quan¬ 
titative  features  differs  to  a  varying  degree  depending  upon  the  approxima¬ 
tions  made  on  aerodynamic  forces  in  the  mathematical  model.  Most  of  the 
comparisons  have  been  made  with  the  linearized  theories  of  flutter,  with  the 
main  emphasis  being  placed  on  the  flutter  boundary  which  delineates,  in  the 
space  of  relevant  paraiaeters,  the  region  wherein  flutter  occurs  from  those 
in  -which  the  shell  is  stable. 

Linearized  models  of  flutter  of  a  circular  cylindrical  shell  have  been 
studied  by  Miles  [8],  Dowell  [9],  Le  Clerc  [10],  and  Voss  [11],  among  others. 
Although  the  linear  theory  is  adequate  for  the  prediction  of  onset  of  flutter, 
it  does  not  provide  any  information  on  the  amplitude  of  the  fluttering 
oscillations.  For  obtaining  such  results,  analysis  of  the  nonlinear  problem 
is  necessary.  Another  important  use  of  the  results  from  nonlinear  analysis 
is  that  on  their  basis  it  can  be  ascertained  as  to  whether  or  not  there  is  a 
possibility  of  the  existence  of  fluttering  motions  at  values  of  the  relevant 
parameters  which  are  lower  than  the  critical  values  predicted  by  linear  theory. 
Finally,  for  circular  cylindrical  shells,  the  linear  theory  fails  to  make  any 
distinction  between  the  circumferentially  traveling  waves  and  standing  waves 
as  modes  of  flutter.  Thus,  although  some  experimental  results  seem  to  suggest 


« 
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that  the  flutter  of  a  circular  cylindrical  shell  is  always  associated  with 
circumferentially  traveling  waves  [5],  the  linear  theory  predicts  that  flut¬ 
tering  motion  in  the  form  of  standing  waves  is  equally  likely  to  occur. 

In  order  to  alleviate  the  aforementioned  limitations  of  linearized 
analysis,  some  analyses  of  nonlinear  flutter  of  a  cylindrical  shell  have 
also  been  conducted  [12,13].  Evensen  and  Olson  [12]  have  essentially  used 
an  assumed  mode  approach  together  with  the  Galerkin's  method  to  analyze 
flutter  of  a  shell  using  shallow-shell  theory  and  first-order  piston  theory 
for  aerodynamic  forces.  Although  their  solution  does  include  both  traveling 
wave  and  standing  wave  responses,  the  accuracy  of  the  results  is  limited  by 
the  fact  that  the  analysis  does  not  include:  (a)  the  nonlinear  terms  in  the 
boundary  conditions  on  axial  stress-resultants  for  the  case  considered  and 
(b)  in-olane  inertia  terms  in  the  equations  of  motion.  Although  one  of  these 
limitations  has  motivated  an  apparently  more  accurate  analysis  [13],  yet  in 
this  solution,  too,  the  nonlinear  boundary  conditions  have  been  satisfied 
only  in  an  averaged  sense,  and,  in  contrast  to  [12],  no  interaction  between 
circumferentially  traveling  waves  is  considered. 

In  what  follows  is  presented  an  analysis  which  takes  into  account  the 
structural  nonlinearities  in  a  consistent  manner.  The  results  obtained  can 
form  a  basis  for  a  computational  algorithm  for  analysis  of  flutter  of  shells 
of  revolution  along  the  lines  of  the  numerical  schemes  for  nonlinear  oscil¬ 
lations  described  in  the  previous  chapters. 


ANALYSIS 


As  for  the  nonlinear  oscillations  problem,  the  analysis  is  most  easily 
carried  out  bv  writing  the  equilibrium  equations  in  the  Koiter's  notation 
[14]  in  the  form 

*~Mli(q,*:q)  +  ~cn(q,fq-  +  -A-^Cq.cq) 

+  P11(q,*q)  +  PnjU.-q)  +  P31(q,-q)  =  o.  O) 

The  vector  q  in  (i)  represents  the  three  displacements  associated  with  the 
shell,  *  is  the  frequencv  of  oscillation,  superposed  dot  denotes  derivative 
with  respect  to  time  nondimensionalized  by  the  frequency  of  oscillation  and 
the  svmmetric  functional  represents  the  inertia  terms.  Further,  the 

terms  P.^Ci  =  1,2,3)  represent  the  first  variation  of  homogeneous  functionals 
of  degree  (i  +  1),  and  they,  of  course,  arise  from  the  structural  stiffness, 
with  the  nonlinearities  being  a  consequence  of  the  quadratic  strain-displace¬ 
ment  relations.  Finally,  the  functionals  and  A^  represents,  respectively, 
the  damping  and  the  nonconservative  forces  on  the  structure.  The  quantity  X 
is  a  "load"  parameter  essentially  representing,  in  case  of  aerodynamic  flutter, 
the  speed  of  the  fluid  flowing  over  the  shell.  Although,  in  general,  the  load 
parameter  can  appear  nonlinear ly  as  well,  and  in  the  other  terms  (particularly 
in  the  damping  term  in  the  equilibrium  equation)  it  is  sufficient  for  describing 
the  basic  methodology  to  assume  the  particular  form  (1). 

It  is  also  noted  here  that  for  shells  of  revolution  each  of  the  com¬ 
ponents  of  q  is  periodic  with  respect  to  the  angular  coordinate:  this  is  a 
result  of  the  invariance  of  the  equilibrium  equations  with  respect  to  angular 
rotation  and  it  leads  to  an  interesting  problem  of  interaction  between  cir¬ 
cumferentially  traveling  waves  in  nonlinear  flutter. 
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Due  to  the  scaling  of  the  time  coordinate  in  (1)  by  the  frequency  of 
oscillation  ,  ve  seek  time  periodic  solutions  of  this  equation  with  a 
period  2~.  The  equation  has  a  trivial  solution  for  all  values  of  a,'.  The 
condition  for  the  existence  of  bifurcation  points  on  the  trivial  solution 
branch  leads  to  the  linear  eigenvalue  problem 


'M,i(y.'q)  +  i^Cin(y,:q)  +  ?  A  Cy.fq)  +  P  (y,  Sq)  =  0, 


(2) 


where  -  is  the  frequency  at  the  onset  of  flutter  and  \  is  the  associated 
o  o 

load  parameter. 

Equation  (2)  is  derived  from  (1)  on  the  assumption  that  the  solution 
at  the  incipient  loss  of  instability  is  of  the  form 

it 


v  e 


so  that  the  functions 

-  -it  -it  -  it 
v  e  ,  v  e  .ye 

are  such  solutions  as  well.  Due  to  rotational  symmetry  the  flutter  mode  v 
is  of  the  form 

y  =  y  e  (3) 

where  y  depends  upon  the  axial  coordinate  only. 

It  is  convenient  to  also  denote  here  the  solution  of  the  problem  adjoint 
to  (2)  by  z,  i.e. ,  z  satisfies 


-woZMn(z,6q)  ”  ^oCjj/z.Sq) 


+  XoA^1(z,6q)  +  P^Cz.iq)  =  0  (4) 

*  * 

where  and  A.^  represent  the  ad  joints  of  the  functionals  and  A^, 
respectively;  thus,  for  any  displacements  q^, 

* 

Cll*ql,q2)  =  Cll(q2’ql)* 

For  the  values  of  (A,w)  near  (*Q»w0)»  the  solution  of  (1)  is  written  as 
q  *  (a  y  e*C  +  c.c.)  +  (8  y  e*C  +  c.c.)  +  W, 


(5) 
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where  w  is  as  yet  an  unknown  function  satisfying  an  orthogonality  condition 


e  1 1  T}  i  ^  * z )  d  t  =  0 , 


e  1C  TuC,:,z)dt  =  0, 


where  1  ^  is  the  first  variation  of  an  arbitrarily  positive  definite  quadratic 

functional.  It  nay  be  noted  here  that  the  choice  of  T  together  with  the 

orthogonality  conditions  (6)  fixes  the  amplitude  a  and  2  uniquely,  i.e. 

2ir 

a  =  ’  e  it:  T11(q,z)dt  (7a) 

f2"‘  -it 

£  =  e  T11(q,z)dt  (7b) 

where  we  have  assumed  that  the  solutions  of  (2,4)  have  been  normalized  so  that 


Tn(y,z)dt  =  1 

'  o 


Tu(y,z)dt  =  1. 


Just  as  the  solution  has  been  decomposed  in  (5),  we  also  decompose  equation 
(1)  into  three  equivalent  ones: 


Ln(W,5q) 


Z7I 

R(q,6q)  -  {e  tT11(y,6q)  e  tR(q,z)dt  +  c.c.) 

n 


it  -  f2n  -it 
-{e  T^Cy,  O  e  R(q , z)dt  +  c.c}  , 


e  lt:R(q,z)dt  =  0, 


|  e-itR(q,z)dt  -  0, 

'  0 


where 


+  Vn(q’£q)  +  pn^q,5q^ 


(12) 
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Once  the 
in  (10) 


solutions  of  (11)  have  been  obtained,  the 
are  calculated  simplv  from  the  relations 


(2)  -  (2) 

W22  =  W11  *  W2 


(o)  _  (o) 


=  w 


11 


other  participating  modes 


which  are  obtained  by  examining  the  nonhotnogeneous  terms  in  equations  for 

(2)  (o)  „ 

w01>  ,  W.70  etc. 

•'ith  (8a)  having  been  solved  in  the  form  of  the  series  (10),  we  can  sub¬ 
stitute  the  result  into  (8b, c)  to  obtain  the  nonlinear  equations  relating 
the  flutter  amplitudes,  frequency  and  the  load  parameters.  The  equations 


turn  out  to  he 

-m(^  -  -)-ji  +  ic(„  -  ^  )a  +  (?  -  X  )aa  +  Y,a2u  +y  =  0,  (13a) 

o  o  o  1  l 

-m(u:"  -  uio2)S  +  ic  (~  -  -Q)t  +  (X  -  +  Y^B26  +  v^Bdt  =  0,  (13b) 

where 

m  =  K^(y,z),  (14a) 

c  =  C11(y,z),  (14b) 

a  =  A11(y,z),  (14c) 

and 

Y1  =  P111^W11(2)’  y,a)  +  Plll(wll(1),y’i)  +  P2ll(y,y’^’  (15a) 

y  2  =  2P2i(y»w12(o))  +  pm(wi2(2),y’i) 

+  pin(wii^0^  *y  >*>  +  2p2ii(y’y^)-  (15b> 


Thus  equations  (13a, b)  are  the  final  results  of  the  analysis.  From  these 
equations  one  can  obtain  both  the  circumferentially  traveling  wave  response 
and  the  standing  wave  response,  as  well  as  the  response  that  arises  due  to 


the  interaction  between  them.  Although  these  equations  are  four  equations 
in  six  variables  —  '• ,  and  the  complex  amplitudes  a  and  £  —  it  can  be 
shown  that  they  admit  solutions  such  that  the  phases  of  the  (complex)  ampli¬ 
tudes  are  arbitrary.  Thus  the  equations  (13)  can  be  used  for  the  calculation 
of  the  lead  parameters  and  frequency  of  harmonic  flutter  as  functions  of 
the  amplitudes  of  the  two  traveling  waves. 

We  end  this  section  bv  emphasising  that  the  results  presented  above  are 
asymptotic  in  nature,  and,  therefore,  they  are  valid  for  sufficiently  small  but 
nonvanishing  amplitudes.  (We  note,  also,  that  the  linear  theory  of  flutter 
is  valid  only  for  the  limiting  case  of  vanishing  amplitude.) 

CONCLUDING  REMARKS 

We  have  presented  an  algorithm  whose  implementation  within  the  frame¬ 
work  of  the  finite  element  methodology  can  provide  a  useful  means  for  ob¬ 
taining  the  effect  of  structural  nonlinearities  on  the  flutter  of  shells  of 


revolution. 
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